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SBLI’s are not trivial in nature 
and are very three 
dimensional flows. 

Physics associated with 
SBLI’s that are often ignored 
in numerical modeling: 

- Heat transfer boundary 
conditions 

- Geometry sensitivities 


Figure from “Computational Fluid Dynamics Investigation 
into the Shock Boundary Layer Interactions in the “Glass 
Inlet” Wind Tunnel” by D. Galbraith, courtesy of M. Galbraith 


- Laminar vs. turbulent flow 
assumptions 
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• Workshop held at the 48 th AIAA Aerospace 
Sciences Meeting. 

- CFD analyses failed to match experimental data. 

• Further CFD analyses performed at the 
University of Cincinnati/NASA Glenn Research 
Center. 

- University of Michigan Glass Tunnel 

• Mach 2.75 freestream 

• 7.75 degree semi-spanning wedge 
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• Focus on the u and v velocity components. 

- Felt that the CFD and post-processing calculations 
from the workshop missed the peak u velocity as well 
as the location of the shock as defined by the v 
velocity profile. 

• Explored alternatives to the previous workshop 
error metric. 
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Geometry and Modeling 

3D overset grid* with 56 million grid points 
divided into 15 zones. 



"Grid based on the one made by Marshall Galbraith 
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• OVERFLOW Version 2.2E 

• Ran on 20 Quad-Core Xeon X5570 (NASA 
Pleiades-Nehalem). 

• Local time-step scaling. 

- CFLMIN=5 

- CFLMAX=20 

• Cases took about 68hrs to converge. 
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CFD Cases 

• Standard 

• Isothermal 

• Modified Geometry 

• Trip 

• Combined 

• TKE 

• MUT 

• Particle Lag* 

• Total Temperature** 

• Perfect vs. Non-Ideal** 
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*Post-processing only 
**Quasi-1 D only 
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• Geometry: As designed (A/A*=3.7062) 

• SST turbulence model 

- Modified SST (SST-GY) 

- BSL 

• All surfaces adiabatic 

• TKE INF =3. 576x1 O' 1 m 2 /s 2 

• R e T,INF = 0-3 

• Constant c p /c v =1 .4 
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• Geometry: As designed (A/A*=3.7062) 

• SST turbulence model 

• Top and bottom walls and wedge isothermal 
(295.7 K), all other surfaces (including bottom 
window) adiabatic. 

• TKE INF =3. 576x1 O' 1 m 2 /s 2 

• R e T,INF = 0-3 

• Constant c p /c v =1 .4 
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Modified Geometry Case 


• Geometry: As currently installed with max 
tolerance (A/A*=3.7847) 

• SST turbulence model 

• All surfaces adiabatic 

• TKE INF =3. 576x1 O' 1 m 2 /s 2 

• R e T,INF = 0-3 

• Constant c p /c v =1 .4 
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• Geometry: As designed (A/A*=3.7062) 

• Laminar from inlet to 67.6mm downstream of the 
throat, SST turbulence model for remaining 
regions. 

- Trip at approximately where Re e =400, based on all- 
laminar case (see next slide) 

• All surfaces adiabatic 

• TKE INF =3. 576x1 O' 1 m 2 /s 2 

• R e T,INF = 0-3 

• Constant c p /c v =1 .4 
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Momentum Thickness Reynolds Number for the Bottom Wall, Center-Span 



X (mm) 
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• Geometry: As currently installed with max error 
(A/A*=3.7847) 

• Laminar from inlet to 67.6mm downstream of the throat, 
SST turbulence model for remaining regions. 

- Ran with SST-GY and BSL in addition to SST. 

• Top and bottom walls and wedge isothermal (295.7 K), 
all other surfaces (including bottom window) adiabatic. 

• TKE INF =3. 576x1 0 1 m 2 /s 2 

• R e T,INF = 0-3 

• Constant c p /c v =1 .4 
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• Geometry: As designed (A/A*=3.7062) 

• SST turbulence model 

• All surfaces adiabatic 

• TKE inf =3.576x 10 3 m 2 /s 2 

• R e T,INF = 0-3 

• Constant c p /c v =1 .4 
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• Geometry: As designed (A/A*=3.7062) 

• SST turbulence model 

• All surfaces adiabatic 

• TKE inf =3.576x 10 3 m 2 /s 2 

• R e T,INF = 3-0 

• Constant c p /c v =1 .4 
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Particle Lag Simulation 



Time constants represent 
50%, 75%, and 1 00% 
total particle relaxation 
time: 

- Short Lag (1 .8 pis) 

- Medium Lag (3.7 ps) 

- Long Lag (5.5 ps) 



Figure 3.1: Measured particle response through an oblique shock. The velocity com- 
ponent normal to the shock, t* n , is normalized by the pre-shock (U ni ) and 
post-shock (u n2 ) velocities and shown as a function of the shock-normal 
direction, n. An exponential fit to the data reveals the particle relaxation 
time, t p — 5.5 /as. 


X = X + UT 



Figure from “Experimental Study of Passive Ramps for Control of 
Shock-Boundary Layer Interactions” by A. Lapsa 
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Total Temperature Sensitivity 


• Discrepancy in total temperature: 

- Workshop: 293 K 

- Experiment: 295.7 ±1 K 

• Using 1 D perfect gas equations: 

- 2.8 m/s (0.47% of a 600 m/s freestream velocity). 

- ±1 K alone is ± 1 m/s (0.1 7% of a 600 m/s freestream 
velocity). 
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Perfect vs. Non-Ideal Air 


• MATLAB code 
developed to perform 
quasi-1 D flow 
calculations for perfect 
and non-ideal air. 

• Very little difference 
between perfect and 
non-ideal air 
calculations 


% Deltas Between Non-Ideal and Perfect Air 
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• 2D Zero Pressure Gradient case from Turbulence Model 
Benchmarking Working Group. 

• Ran with SST, SST-GY, BSL, and K-Omega 
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Diagram from http://turbmodels.larc.nasa.gov/FlatPlate/ 
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Normalized Turbulent Shear Stress at Re=4.2E6 




Solutions agree well with each other (except K-Omega) 
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Resolution Study 




Inlet: Tt=295.7 K, Pt=98000 Pa 


Throat 



Coarse 

Medium 

Fine 

Tt (K) 

296.530 

295.909 

295.701 

Pt (Pa) 

98761 .8 

98199.4 

98009.4 

M 

0.95217 

0.94737 

0.94703 


Upstream 



Coarse 

Medium 

Fine 

Tt (K) 

296.504 

295.960 

295.704 

Pt (Pa) 

98134.4 

98103.0 

97996.9 

M 

2.73365 

2.73637 

2.74482 


Downstream 



Coarse 

Medium 

Fine 

Tt (K) 

296.173 

296.073 

295.742 

Pt (Pa) 

96083.7 

96459.3 

96247.1 

M 

2.49111 

2.48308 

2.47317 


Mass flow conserved within < 0.5% 
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Results 


www.nasa.gov 23 


National Aeronautics and Space Administration 


w 

Cincinnati 


UNIVERSITY OF 



Data Comparison Plane 
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Standard Case CFD 


U/UInf 
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Standard Case CFD, with Short Lag 


U/UInf 
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Standard Case CFD, with Short Lag 
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Standard Case CFD, with Medium Lag 


U/UInf 
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Standard Case CFD, with Medium Lag 
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Standard Case CFD, with Long Lag 


U/UInf 
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Standard Case CFD, with Long Lag 
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Isothermal Case 
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Static Temperature at X=18.191mm 
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Positive U-Velocity Difference 
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X (mm) 


Negative U-Velocity Difference 
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Isothermal case shifts interaction region slightly upstream 

Difference=lsothermal-Standard 
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Case 
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CFD turbulent kinetic energy lower in the freestream 
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TKE and MUT Cases 
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Increase in freestream TKE results in better agreement with 

freestream experimental data 
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Turbulence Model Effects 
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Metrics 
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Note all prior workshop 
CFD analyses utilized a 
total temperature of 293 K 
while the new CFD 
analyses utilized a total 
temperature of 295.7 K. 


U Error 

V Error 

0.02373 
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F 
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Combined (Medium Lag) 
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M 
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1 
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Point Comparison Location 


A: Taken at center- B: Point of most upstream 
height, center-span U max within experimental 

data plane 
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Point A 

Point B 

Case 

U (m/s) 

AU (m/s) 

U (m/s) 

AU (m/s) 

Standard 

594.600 

0.000 

587.042 

0.000 

Isothermal 

594.454 

-0.146 

586.377 

-0.665 

Modified Geometry 

596.567 

1.967 

586.413 

-0.629 

Trip 

595.186 

0.586 

587.737 

0.695 

Combined 

596.980 

2.380 

587.038 

-0.004 

Experiment 

- 

- 

599.330 

1 2.288 


A=Case-Standard 
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• CFD analyses were performed and generally 
under predicted the freestream velocities but 
with improvements. 

- Improved modeling. 

- The flow was shown to be most likely transitional 
downstream of the throat. 

- SST likely has corner separation too large, which was 
reduced with SST-GY and BSL. 
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• A fraction of the measured PIV lag was used 
with a simple model to modify the CFD solutions. 

- Showed improved comparisons to the experimental 
data. 

- Future comparisons should have the CFD results 
augmented in a post-processing step to calculate 
particle velocities. 

• New complimentary metrics: 

- max u velocity 

- shock angle 
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• Sensitivities to address: 

- Additional geometric parameters 

- Turbulence model and parameters 

- Heat transfer boundary conditions 

• Conjugate heat transfer 

- Boundary-layer transition/trip location and model 
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• The simplified PIV model should be improved 
on. 

- Calculating the particle lag based on the forces 
exerted on the individual particles by the air (including 
particle size distribution). 

- Obtaining flow field snapshots at two instances in 
time. 

• Snapshots would then be processed using the same PIV 
post-processing algorithm used with the experimental data. 
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Throat (in) 

Test Section (in) 

A/A* 

As Designed 

2.25 x 0.742 

2.25 x 2.75 

3.7062 

As Installed 

2.25 x 0.725 

2.25 x 2.72 

3.7847 


Error of “As Installed” Measurement: +/- 0.005 in 
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Throat Modification 



Added an additional 50 grid points to base grid to define the nozzle 

contour 
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• Raised bottom wall. 

• Modified nozzle 
contour. 


y new 



,new 


y throat. 


X x -X 
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1 . Throat 

2. Trip Location 

3. Start of Straight Section 

4. Wedge Leading Edge 

5 . Wedge Trailing Edge 


Tt(K) 
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Total Temperature at Bottom Wall, Center-Span 
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